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PREFACE 

This  Memorandum,  prepared  for  the  Advanced  Research  Projects 
Agency;  is  part  of  a  study  of  the  effect  of  atmospheric  turbulence 
on  optical  and  infrared  reconnaissance  and  guidance  equipment. 

A  quantitative  understanding  of  the  manner  in  which  an  initially 
coherent  beam  of  finite  cross  section  propagates  is  required  for  the 
prediction  of  the  performance  of  various  devices  employing  lasers 
for  target  acquisition  or  guidance  in  tactical  missions,  optical 
communication  systems,  and  other  devices.  This  Memorandum  calculates 
the  mean  intensity  distribution  for  an  arbitrary  amplitude  and  phase 
distribution  in  a  finite  aperture  in  both  the  near  and  far  field  and 
examines  in  detail  the  case  of  a  uniform  distribution  across  a  circu¬ 
lar  aperture. 

These  results  should  be  of  use  to  those  interested  in  propa¬ 
gation  theory  and  its  applications  to  laser  range  finders,  laser  \ 

line  scanners,  communication  systems,  and  various  guidance  and  other 
systems  employing  an  illuminating  beam. 
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SUMMARY 

The  £irst  part  of  this  Memorandum  is  devoted  to  extending  the 
Huygens -Fresnel  principle  to  a  medium  which  exhibits  a  spatial  (but 
not  temporal)  variation  in  index  of  refraction.  With  the  proof  of  a 
reciprocity  theorem  for  a  monochromatic  disturbance  in  a  weakly  inhomo¬ 
geneous  medium,  it  is  shown  that  the  secondary  wavefront  will  be  deter¬ 
mined  by  the  envelope  of  spherical  wavelets  from  the  primary  wavefront, 
as  in  the  vacuum  problem,  but  that  each  wavelet  is  now  determined  by 
the  propagation  of  a  spherical  wave  in  the  refractive  medium. 

In  the  second  part,  the  above  development  is  applied  to  the  case 
in  which  the  index  of  refraction  is  a  random  variable;  a  further 
application  of  the  reciprocity  theorem  results  in  a  formula  for  the 
mean  intensity  distribution  from  a  finite  aperture  in  terms  of  the 
complex  disturbance  in  the  aperture  and  the  modulation  transfer 
function  (MTF)  for  a  spherical  wave  in  the  medium.  The  results  are 
applicable  in  both  the  Fresnel  and  Fraunhofer  regions  of  the  aperture. 
Using  a  Kolmogorov  spectrum  for  the  index  of  refraction  fluctuations 
and  a  second -order  expression  for  the  MTF,  the  formula  is  used  to 
calculate  the  mean  intensity  distribution  for  a  plane  wave  diffracting 
from  a  circular  aperture  and  to  give  approximate  expressions  for  the 
beam  spreading  at  various  ranges.  It  is  argued  that  the  spherical- 
wave  MTF  is  the  basic  quantity  to  be  measured  for  computing  the 
atmospheric  degradation  of  an  intensity  pattern  from  an  arbitrary 
disturbance  in  a  finite  aperture. 
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I.  INTRODUCTION 

The  problem  of  computing  the  mean  diffraction  pattern  from  an 
optical  wave  in  a  finite  aperture  in  the  presence  of  a  turbulent 
atmosphere  is  both  of  considerable  Importance  for  a  number  of  prac¬ 
tical  applications  (e.g.,  laser  radar  and  laser  line  scanners)  and 
a  useful  tool  for  evaluating  our  understanding  of  turbulence 
theory. 

To  solve  the  problem,  it  is  first  shown  in  Section  II  that 
the  Huygens -Fresnel  principle  may  be  extended  to  a  refractive  medium, 
i.e.,  the  field  due  to  some  disturbance  specified  over  an  aperture 
can  be  computed  by  superimposing  spherical  wavelets  radiating  from 
all  elements  of  the  aperture.  The  principle  follows  directly  from 
Green's  theorem  and  the  Klrchhoff  approximation  applied  to  the  scalar 
wave  equation,  together  with  a  field  reciprocity  theorem  (proven  in 
the  Appendix)  between  an  observation  point  and  a  source  point  of 
spherical  waves  in  an  Inhomogeneous  medium.  To  apply  the  extended 
principle  to  the  atmosphere,  we  assume  that  the  refractive  index 
fluctuation  is  a  random  function  of  coordinates  and  does  not  depend 
on  time,  i.e.,  the  refractive  index  n  does  not  change  appreciably 
over  periods  of  the  order  At  ■  >/c  (At  ■  3  x  10  ^  sec  for  a  wave¬ 
length  \  ■  2rr/k  of  lp) ,  where  c  is  the  vacuum  speed  of  light.  The 
time  changes  in  n  are  considered  as  changes  in  the  different  realiza¬ 
tions  of  the  random  field  n(r).  On  the  basis  of  this  principle,  one 
can  separate  the  geometry  of  the  problem,  i.e.,  the  disturbance  in 
the  aperture,  from  the  propagation  problem,  which  is  determined  by 
the  manner  in  which  a  spherical  wave  propagates  through  the  medium. 
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The  mean  intensity,  calculated  in  Section  111,  is  found  by  first 
computing  the  intensity  at  a  point  from  an  arbitrary  pair  of  elements 
in  the  aperture.  Applying  the  reciprocity  relationship  and  averaging 
reveal  that  the  above  quantity  is  essentially  the  modulation  transfer 
function  (MTF)  for  a  spherical  wave  in  the  medium.  The  integration 
over  the  aperture  is  performed  as  a  final  step,  which  results  in  a  sim¬ 
ple  formula  for  the  mean  intensity  pattern  in  both  the  Fresnel  and 
Fraunhofer  regions  for  an  arbitrary  complex  disturbance  in  the  aperture. 
The  properties  of  the  medium  (e.g.,  the  turbulence  parameters)  appear 
only  in  the  MTF.  The  MTF  describes  the  reduction  in  lateral  coherence 
between  different  elements  of  the  aperture,  effectively  transforming 
the  aperture  into  an  equivalent  partially  coherent  radiator,  with  the 
degree  of  coherence  decreasing  with  increasing  distance  from  the 
aperture.  As  a  practical  example,  the  mean  intensity  patterns  are 
computed  over  horizontal  ranges  of  interest  (of  the  order  of  kilometers) 
for  both  a  2 -cm  and  6-cm  aperture  using  a  range  of  turbulence  parameters 
characteristic  of  paths  of  the  order  of  a  few  meters  above  the  ground. 
Over  certain  limited  path  lengths,  approximate  formulas  exist  for  the 
MTF  for  a  Kolmogorov  spectrum.  Introducing  these  expressions  into 
the  above-mentioned  formula  yields  approximate  expressions  for  the 
beam  spreading  at  these  ranges.  It  is  shown  that  at  distances  suf¬ 
ficiently  large  so  that  the  coherence  length  of  the  MTF  is  less  than 
the  size  of  the  smallest  inhomogeneities ,  the  Intensity  pattern  is 
dominated  by  multiple  scattering  from  these  smallest  scatterers  and 
is  contained  in  an  angle  which  is  independent  of  both  wavelength  and 
aperture  size. 
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It  Is  also  shown  that  (1)  for  homogeneous  isotropic  turbulence 
an  optical  measurement  of  the  spherical -wave  MTF  can  be  inverted  to 
give  the  turbulence  spectrum,  and  (2)  the  intensity  distribution  can 
be  Inverted  to  give  the  spherical-wave  MTF,  but  only  for  spatial 
separations  smaller  than  the  diameter  of  the  transmitting  aperture. 

Finally,  it  is  argued  that  the  spherical-wave  MTF  is  the  basic 
quantity  to  be  measured  for  computing  the  degradation  of  an  intensity 
pattern  from  an  arbitrary  disturbance  in  a  finite  aperture. 


II.  AN  EXTENSION  OF  THE  HUYGENS -FRESNEL  PRINCIPLE 
TO  A  REFRACTIVE  MEDIUM 

The  problem  considered  in  this  section  is  the  calculation  of 
the  complex  field  at  an  arbitrary  point  in  a  turbulent  medium,  given 
a  complex  monochromatic  disturbance  specified  over  a  finite  surface 
(i.e.,  the  aperture). 

We  assume  a  scalar  wave  equation  for  the  propagation  of  an 
electromagnetic  disturbance  in  the  optical  wavelength  region  through 
a  nonabsorbing  refracting  medium: 

(V2  +  k2n2)  U  -  0  (1) 

where  U  is  one  component  of  the  electric  or  magnetic  field  vector 
ard  n(r)  ■  1  +  n^Cr)  is  the  refractive  index  of  the  medium  at  the 
point  :r,  where  (n)  ■  1  and  n^  is  the  small  fluctuating  part  of  the 
refractive  index  field.  (Angular  brackets  are  used  to  denote  ensemble 
averages .) 

The  method  to  be  used  is  analogous  to  the  integral  theorem  of 
Helmholtz  and  Kirchhoff^  used  in  vacuum  diffraction  theory,  and  a 
similar  notation  wil]  be  employed. 

Let  the  turbulent  medium  occupy  a  volume  V  and  be  bounded  by  a 
closed  surface  S,  and  assume  that  U  possesses  continuous  first-  and 
second-order  partial  derivatives  within  and  on  S.  Within  V,  let  U/ 
be  the  field  at  £  due  to  a  point  source  at  P.  Then  U'  satisfies  the 


equation 
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where  6  is  the  Dirac  delta  function.  Multiplying  Eq.  (1)  by  U7, 

Eq.  (2)  by  U,  subtracting,  integrating  over  V,  and  applying  Green's 
theorem  yields 


U(P) 


[UVU7  -  u'vul  •  dA 


(3) 


where  dA  denotes  the  surface  element  with  normal  directed  into  V. 

To  proceed  further,  we  apply  the  Kirchhoff  approximation.  A 

monochromatic  disturbance  is  assumed  to  propagate  through  an  opening 

in  a  plane  opaque  screen  which  is  large  compared  to  the  wavelength 

but  small  compared  to  the  distance  between  P  and  the  screen.  The 

disturbance  at  P  is  computed  by  taking  Kirchhoff 's  integral  over  a 

surface  S  formed  by  (1)  the  opening  A,  (2)  a  portion  B  of  the  non- 

illuminated  side  of  the  screen,  and  (3)  a  portion  C  of  a  large  sphere 

of  radius  R,  centered  at  P,  which  together  with  A  and  B  forms  a  closed 

surface.  The  contribution  to  the  field  at  P  from  C  can  be  argued  to 

vanish  by  assuming  that  the  radiation  field  does  not  exist  for  all  times, 

but  begins  to  radiate  at  t  ■  t  (which  actually  implies  an  unimportant 

departure  from  monochromacy) .  Then  if  for  all  times  of  interest  R  is 

taken  to  be  >  c(t  -  tQ) ,  the  contributions  from  the  surface  C  could 

* 

not  have  reached  the  point  P,  and  the  integral  will  vanish.  Finally, 

invoking  the  Kirchhoff  approximation,  we  assume  that,  everywhere  on 

3U  3U 

B,  U  and  —  will  be  approximately  zero  and  that,  on  A,  U  and  will 

not  differ  appreciably  from  the  values  obtained  in  the  absence  of 

*  7l) 

Born  '  has  shown  that  this  assumption  is  not  essential,  but  it 
shortens  the  discussion. 


* 
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the  screen.  Hence  the  surface  Integral  In  Eq.  (3)  reduces  to  an 
Integral  over  the  aperture  area,  A,  and  from  this  assumption  one  can 
develop  an  extension  to  the  Huygens-Fresnel  principle  for  vacuum  prop¬ 
agation  to  propagation  In  a  weakly  Inhomogeneous  medium.  The  princi¬ 
ple  as  applied  to  propagation  in  a  vacuum  states  that  every  point  of 
a  wavefront  may  be  considered  as  a  center  of  a  secondary  disturbance 
which  gives  rise  to  spherical  wavelets,  and  the  wavefront  at  a  later 
Instant  may  be  regarded  as  the  envelope  of  these  wavelets.  It  will 
be  shown  that,  for  a  refractive  medium,  the  "extended"  principle  Is 
that  the  secondary  wavefront  will  again  be  determined  by  the  envelope 
of  spherical  wavelets  from  the  primary  wavefront,  but  each  wavelet 
will  now  be  determined  by  the  propagation  of  a  spherical  wave  through 
the  refractive  medium.  In  either  case,  the  theorem  Is  an  approxima¬ 
tion  valid  only  when  the  scattering  angles  are  sufficiently  small. 

We  first  assume  the  aperture  to  lie  In  a  plane  normal  to  the 

z  direction,  and  we  note  that  when  the  radius  of  curvature  at  each 

3U 

point  of  the  aperture  Is  large  compared  to  the  wavelength,  ^  «w 
lkU(e  •  e  ),  where  e  Is  the  unit  vector  normal  to  the  wavefront  In 

™ 2  "TV  •  |  i  I 

Iks  +  t 

the  aperture.  Then,  writing  U7  ■  6|  1  1  ’■■  ,  where  s  ■  |r  -  p|  is  the 

geometric  distance  between  the  point  F  and  the  elemental  area  In  the 
aperture,  and  assuming  ks  »  1  yield 


ik  r  elks  +  ♦  r 

-  ft  JA  "(*>  — —  [• 


e  •  e  -  e 
—  — z  — n 


'  +  lk  •  *J  d2t  (4> 


where  d  r.  Is  an  element  of  area  at  the  point  r  in  the  aperture  and 
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e  is  a  unit  vector  in  the  direction  from  P  to  r.. 

8  iks  +  * 

For  r,  ^  P ,  substituting  U7  ■  - - into  Eq.  (2)  yields  the 

differential  equation  for  *: 

V2*  +  (V*)2  +  2V*  •  V*q  +  k2(2nL  +  n2)  «  0  (5) 


where 


*  iks 
o  m  e 

s 

Expanding  *  *  ,  where  ^  ~  n*,  in  Eq.  (5)  generates  the 

(3) 

hierarchy  of  Rytov  solutions  to  the  field  at  r  from  a  unit  spherical 

wave  emitted  at  P  in  an  inhomogeneous  medium.  In  the.  geometric  optics 

regime,  *  ~  ik  J  n^ds,  •  e^l  as  In^^  •  e^)  | ,  and  in  the  atmosphere 

at  optical  frequencies,  |n^|  ~  10  **.  In  general,  due  to  the  smallness 

of  both  the  scattering  "potential"  (proportional  to  n^)  and  the 

ratio  of  the  wavelength  to  the  size  of  the  smallest  inhomogeneity  (of 

dimension  Lq) ,  from  a  photon  picture  the  scattered  wave  from  the  point 

2 

source  will  vary  slowly  over  a  wavelength  for  as  (kXQ)  scatterings. 

The  scattering  length  can  be  estimated  from  the  decay  rate  of  the 

average  field  for  a  single-scale  model  to  be  as  ((n^)  k  ,  which 

implies  that  the  V*/k  term  can  be  neglected  compared  to  unity  for  all 

2  6 

distances  of  interest  (i.e.,  I  /( n, )  as  10  km). 

o  l 

It  might  be  noted  that  at  a  range  where  |v*/k|  ~  1  one  would 
measure  a  local  beam  divergence  of  the  order  of  a  radian  due  to  the 
scattering,  which  has  never  been  observed. 

Then,  noting  that  in  our  development  the  area  A  could  have  been 
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replaced  by  any  other  open  surface  whose  rim  coincides  with  the  edge 
of  the  aperture,  we  choose  instead  of  A  a  portion  of  the  incident 
wavefront  which  approximately  fills  the  aperture  and  neglect  any  small 
errors  near  the  rim.  Over  the  new  region  of  integration,  A',  the  de¬ 
rivative  of  U  in  the  direction  normal  to  the  surface  is  =•  ikU,  and 
defining  x  “  tt  -  cos  *  (e^  •  e^),  we  obtain 

U(P)  =  [  K(x)  U'(r,P)  U(A')  dA'  (6) 

JA/ 

where  K(x)  “  (-i/2\)(l  +  cos  x)  • 

In  the  Appendix  we  show  that  U'(r,P)  =  U  7(P  ,r} ,  i.e.,  reciprocity 

holds  in  that  the  field  and  source  points  may  be  interchanged  in 

Eq.  (2).  Equation  (6)  is  thus  the  extension  of  the  Huygens  principle, 

where  the  contribution  from  the  element  dA '  of  the  wavefront  to  the 

iks  +  f 

field  at  P  is  K(x)  U(A')  - -  dA'. 

s 

When  the  distance  from  the  observation  point  to  the  aperture  is 

large  compared  with  the  linear  dimensions  of  the  aperture  and  the 

wavefront  does  not  vary  appreciably  over  the  aperture,  K(x)  may  be 

replaced  by  K(XQ) >  where  xq  is  the  angle  between  the  normal  to  the 

screen  and  a  line  from  the  center  of  the  aperture  to  the  observation 

point.  In  the  same  approximation,  the  factor  l/s  may  be  replaced  by 

the  reciprocal  of  the  distance  from  the  center  of  the  screen  to  the 

point  P,  1/s  .  Then  Eq.  (6)  reduces  to 
o 

K(X0>  P  iks  +  t  ,  .  2 

-  e  T  U(r)  d  r 

8o  J  "  " 


U(P)  = 


(7) 


Equation  (6)  or  (7)  reveals  that  If  we  know  hew  a  spherical  wave  prop¬ 


agates  In  a  given  medium,  we 
disturbance  In  an  aperture, 
ample  of  the  utility  of  this 


can  determine  the  response  to  an  arbitrary 
In  the  next  section  we  will  give  an  ex¬ 
formulation  by  applying  it  to  a  problem 


of  practical  Importance. 


-10- 


III.  THE  MEAN  INTENSITY  FROM  A  DIFFRACTING  APERTURE 


From  Eq.  (7),  the  intensity  at  the  point  P  is 


I(P) 


x 


♦  (£.)  +  t*(E2>  2 

e  U(£l)  U*^)  d^ 


(8) 


where  s^,  s^  are  the  geometric  distances  between  the  point  P  and  the 
points  r^,  £2  in  the  aperture,  respectively,  and  ir(r^),  are  the 

perturbations  in  the  field  at  P  due  to  unit  spherical  waves  emitted 
at  £2*  For  the  case  in  which  the  index  of  refraction  fluctuations, 
n^,  is  a  random  variable,  the  mean  intensity  from  a  finite  diffracting 
aperture  can  be  computed  if  ^exp  Ct(r^)  +  is  known  for  all 

£^,  £2 •  Similarly,  higher  moments  of  I(P)  can  be  computed  from  a 
knowledge  of  terms  of  the  form  ^exp  Ct(£j_)  +  +  t (r^)  +  f*(r^)]^, 

etc.  For  the  mean  intensity  problem,  it  follows  directly  from  the  re¬ 
ciprocity  relationship  (in  the  Appendix)  that  exp  Ct(r^)  +  ^(r^)]  is 
the  field  at  £^  multiplied  by  the  complex  conjugate  of  the  field  at 
£2  (normalized  to  the  vacuum  field)  due  to  a  spherical  wave  at  £.  In 
particular,  the  quantity 

ik^  -  s2) 

—  («» »<*!>  +  ♦*<%>') 

is  equal  to  the  cross-correlation  of  the  complex  fields  at 


-li¬ 


the  points  r^,  due  to  a  unit  point  source  at  P_.  The  function 
^exp  [i|r(r^)  +  by  definition,  the  MTF  for  a  spherical 

wave,  Ms(rlf  £2>  z) . 

When  one  observes  the  intensity  as  a  function  of  time,  one  sees 

the  effects  of  beam  "wander"  and  "breakup,"  both  of  which  contribute 

to  the  intensity  scintillation  at  a  point.  If  the  average  is  taken 

over  times  long  compared  with  the  periods  over  which  varies,  this 

mean  intensity  distribution  would  include  all  of  the  above  effects. 

If  averages  over  an  ensemble  corresponding  to  a  shorter  time  period 

could  be  obtained,  a  time  history  of  the  intensity  distribution 

corresponding  to  observations  over  this  shorter  time  interval  could 

be  obtained.  For  plane  waves,  the  long-term  MTF,  correct  to  second 

(4) 

order  in  n1 ,  is  given  by 


M  (p,  z)  =  exp 
P  (Zc 


Jo  J0(Kp)  VK)  KdK 


Jo  *n(K)  KdK 


(9) 


where  p  =  |r^  -  I  »  $n00  is  the  spectrum  of  the  index  of  refraction 
fluctuations,  and 


z 

c 


(10) 


is  the  propagation  distance  in  which  the  mean  field  of  a  plane  wave 
or  spherical  wave  is  reduced  by  e  ^  from  its  vacuum  value.  The 
quantity  z is  also  the  distance  over  which  the  perturbations  in 
the  field  due  to  the  medium  become  comparable  to  the  field  in  the  ab¬ 
sence  of  the  medium,  and  can  be  computed  by  setting  1/2^|i|TjJ^  =  1. 

It  should  be  noted  that,  independent  of  the  form  of  the  spectrum, 


zc°-  \  Consistent  with  the  approximation  that  the  angles  x  of 
interest  remain  small  is  the  assumption  that  the  points  r^,  r^ 
may  be  considered  to  be  lying  on  a  large  sphere  centered  at  P.  In 
this  case,  the  MTF  given  by  Eq.  (9)  is  modified  for  spherical  waves 
by  substituting  for  the  Bessel  function  JQ(Kp)  the  quantity 

Jo  J0(Kpu) 

Using  Eq.  (9)  for  the  case  of  homogeneous  isotropic  turbulence, 
one  can  determine  the  turbulence  spectrum  from  optical  measurements 
of  the  spherical  wave  MIF  (e.g.,  by  using  a  sufficiently  small 
aperture) . 

First,  the  distance  Z£  is  determined  by  considering  a  sufficiently 

short  path  length,  z,  and  measuring  the  limiting  value  of  Mg(p,  z)  for 

spatial  separations  much  greater  than  the  largest  turbulence  scale,  L  . 

"2z/zc  ° 

From  Eq.  (9),  this  limiting  value  is  e 

Then,  from  the  relation  between  the  spherical  and  plane  RTFs, 
one  can  construct  the  quantity  ln[Mp(p,  z)],  which  is  given  by 

ln[Mp(p,  z)]  -  ^  [p  In  Ms(p,  z)]  (11) 

With  z£  and  the  logarithm  of  the  plane -wave  MTF  known,  one  can 
invert  Eq.  (9)^  to  yield 


»„(K> 


2  2 
2tt  k  z 


od  2 

J  JQ(Kp)  [l  +  57  In  Mp(p,  z)]  pdp  (12) 


for  the  turbulence  spectrum. 

(A) 

Yura  and  Lutomirski 


have  used  the  spectral  density 
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»n(K)  - 


0.033  C2  g'CKl0/5*92) 

(/♦l'2)11'6 

o 


(13) 


an  extrapolation  of  the  Kolmogorov  spectrum.  In  Gq.  (13),  Is  the 

Index  structure  constant  and  l  (■  2rr i  )  and  L  (■  2ttL  )  are  the  Inner 

oo  o  o 

and  outer  scales  of  turbulence,  respectively.  In  Figs,  la  and  lb  we 

plot  zc  versus  X  for  three  values  of  roughly  corresponding  to  weak 

(3  x  10"16  cm‘2/3),  medium  (3  *  Iff15  cnf2/3),  and  strong  (3  *  10-1A 
-2/3 

cm  )  turbulence.  For  horizontal  propagation  near  the  ground,  we 

have  used  the  nominal  values  of  i  1  0.1  cm  and  t  “  10  cm  and  100  cm. 

o  o 

For  jfcQ  «  f.Q,  one  can  approximate  the  Integral  in  Eq.  (10)  to  yield 


/a  /  u2r2  -  5 / 3 N -1 
z  ^  (0*4  k  C  I  ) 
c  no 


(14) 


The  coherence  length  z  ,  and  hence  the  MTF,  will  depend  strongly  on 

c 

the  outer  scale  of  turbulence. 

Returning  to  Eq.  (8),  we  thus  obtain 


<‘>  ‘  *2  If 


Ik (  s  ^  -  s2) 


Mg(r.1»  £2 *  d2~l  d2-2 


2  2 

where  A  ■  [(1  +  cos  x  )/2s  X]  and  M  Is  the  MTF  for  spherical  waves. 

O  O  8 


For  the  case  of  a  plane  wave  In  a  circular  aperture  of  diameter  D, 
we  take  the  origin  of  the  coordinate  system  at  the  center  of  the 
aperture,  let 
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A 

1 


U<tl,2)  ■  k ,2 1  S  D/2 
"  °>  |rlj2l  >  °/2 


(16) 


and  integrate  over  the  entire  plane  of  the  aperture.  If  we  define 
the  vector  £  as  the  normal  from  the  (z)  axis  of  symmetry  to  the 
observation  point  P,  then 


1,2 


'  +  •  El,2>' 


2 

and,  in  the  small-angle  approximation,  «  (j^-)  and 


-2 2  *  (% 

S1  ’  s2  “ 

Changing  variables  to  £  *  X 

of  homogeneous  isotropic  turbulence 


»  ,  2  2 
-  r2)  +  -  r2 

2z 


(17) 


=  +  £2)  ,  then  for  the  case 


=  (~^)  J  d2£.Mg(p,  z)  e  (ik//z)£.,£.  j*  U(r  +  j£) 


x  U-(r  -  \  s)  e(lk/z)£-£  d2r 


(18) 


For  our  problem  U  is  real,  and  inspection  of  Eq .  (18)  reveals  that 

the  inner  integral  is  the  integration  of  the  function  e(*k/z)fi.*I. 

over  the  area  of  overlap  of  two  circles  each  of  diameter  D,  with 

centers  located  a  distance  p  apart.  The  integration  is  straightforward 

2 

and  yields  for  the  inner  integral  D  r\,(x),  where 

P 


* 
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.C08_1(X) 

r»w  '  J.  { 


sin  [2B.x(cos  6  -  x)  ] 
(2px  cos  6) 


{l  -  cos  -I2j3x(  cos  6  -  x)]} 
(20x  cos  0)2 


} 


d6 


x  s  1 


-  0,  X  >  1 


(19) 


where  x  ”  p/D  and  0  ■  kD  /4z.  Then,  using  polar  coordinates  for  the 
£  Integration,  performing  the  angular  Integral,  and  changing  variables 
from  p  to  x  ■  p/D  yield 


<I>  ■  ”  @2  J  x  Jo(2«x)  Mg(Dx,z)  r„(x)  dx 


(20) 


kDp  kD 

where  J  Is  the  zero-order  Bessel  function  and  a  ■  ’r_K  "  rr  tan  0, 
o  l  z  z 

where  0  («  XQ)  is  the  angle  which  the  direction  to  P  makes  with  the 
central  direction.  Normalizing  the  Intensity  at  a  particular  range 
(z)  to  the  on-axls  value  ( a  “  0)  In  the  absence  of  turbulence  (M  ■  1), 
we  obtain 


x  Jp(2crx)  Ms(Dx,z)  fg(x)  dx 

^  Si  x  rR(x)  dx 


(21) 


2 

In  the  Fraunhofer  region  of  the  aperture  0  -  0,  and  D  ro(x)  re¬ 


duces  to  the  overlap  area  of  the  two  circles.  From  Eq.  (19) 

cos  *(x)  -  x 


r0(x)  -  f 


_1/  v  ^ 


(22) 
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and  substituting  In  Eq.  (20)  with  M  -  1  and  integrating  by  parts 
yields  the  vacuum  Airy  pattern: 


!„(“>  -  7  ** 


2^(0) 


For  arbitrary  3,  Eq.  (21)  gives  the  normalized  mean  Intensity  distri¬ 
bution  in  the  Fresnel  region  as  well  as  in  the  Fraunhofer  region  of 
the  aperture.  Numerical  calculations  show  that  Eq.  (22)  is  a  good 
approximation  to  T^(x)  for  3  ~2.  (See  Fig.  2.) 

We  have  numerically  integrated  Eq.  (21)  for  a  2-cm  and  6-cm  ap¬ 
erture,  using  the  spectrum  of  Eq.  (13)  for  propagation  paths  of  1  and 
5  km  for  X  -  0.6328y  and  5  and  10  km  for  X  -  10. 6y.  The  results  are 
shown  in  Figs.  3  through  6  for  turbulence  parameters  characteristic 
of  the  strengths  and  scales  found  from  zero  to  several  meters  above 
the  ground.  The  dashed  curve  on  each  graph  is  the  vacuum  intensity 
pattern  at  the  range  indicated,  normalized  to  its  value  at  0-0.  For 
a  given  range  the  6-cm  aperture  has  less  effective  coherence  than  the 
2-cm  aperture,  which  results  in  a  correspondingly  greater  degradation 
of  the  distribution  in  the  absence  of  turbulence.  Further,  for  given 
range,  aperture  size,  and  turbulence  conditions,  the  longer  wavelength 
10. 6y  pattern  is  the  more  coherent  due  to  the  larger  coherence 
length,  zc< 

In  Fig.  7a-c  we  compare  the  mean  intensity  distributions  for 
the  above  wavelengths  and  aperture  diameters.  As  explained  in  the 
next  section,  the  angle  for  which  the  intensity  is  reduced  to  one- 
half  its  axial  value  becomes  independent  of  both  X  and  D  at  large 


ranges . 


* 


Fig.  2  —  The  overlap  integral  defined  by  Eq.  (19) 


A  *1.06/4,  0  *2  cm 
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IV .  APPROXIMATE  FORMULAS 


The  integrals  involved  in  Eq.  (20)  for  computing  the  mean  pattern 
can  be  economically  performed  on  a  computer.  However,  for  the  purpose 
of  relating  our  work  to  others  and  to  show  precisely  what  conditions 
are  required  for  their  use,  formulas  for  the  beam  spread  over  certain 
limited  ranges  are  obtained. 

If  p  is  the  transverse  distance  for  which  M  (p  ,z)  =  e  ^ .  then 
for  the  modified  Kolmogorov  spectrum  of  Eq .  (13),  approximate  formu¬ 
las  exist  for  the  MTF  when  p  «  l  and  i  «  p  «  L  .  (It  should  be 

0  0  o  o  o 

observed  that  the  condition  |vty/k|  «  1  implies  a  lower  limit  of 
Pq  ~  X . )  Introducing  these  MTF  a  into  Eq .  (20)  results  in  approximate 
formulas  for  the  mean  intensity,  provided  that  the  correlation  length 
Pq  in  each  case  is  small  compared  with  the  size  of  the  aperture. 

CASE  A:  p  «  a 
o  o 

When  sufficient  scattering  has  occurred  so  that  the  mutual 
coherence  between  the  field  from  different  elements  of  the  aperture 
is  small  compared  with  the  smallest  scale  of  turbulence,  SLq1  the 
Bessel  function  in  Eq .  (9)  can  be  expanded  in  powers  of  (p llQ)>  and 
the  MTF  for  a  spherical  wave  can  be  approximated  by 


Mg(p,  z)  =  exp 


I  zo2  -To  VK> R3dK  I 
I  6z°  Jo  VK>  KdK  I 


=  exp 


,  2  2 

k  zp  q 


} 


(23) 


where  q  =  tt^/3  ^  $  (K)  K^dK. 
•J  u  n 


Equation  (22)  will  be  valid  for  all  p's 
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of  interest  if  and  only  if 


z  »  z,  = 


2  2  -1 

(kV  q)  1 


For  the  spectrum  given  by  Eq.  (13),  q  can  be  estimated  by  computing 


- (K£  / 5.92) 
„3  o 


I*  K  e  dK  =  1  (1^2)  r 

^0  (K2  +  i‘2)ll/6  2  '  ^o  ■’( 


0  [y  +  (£  /5 .92  £  )*] 
o  o  J 


2.11/6 


■HW  ra/«^)' 


where  T  is  the  gamma  function  and,  for  &  «  L  , 

o  o 


q  0.5  C2 

no 


2  2  2 

Defining  y  =  k  D  zq,  it  follows  from  Eq .  (24)  that  for  D  SL 

2 

y  »  1.  When  M  (Dx ,  z)  is  substituted  into  Eq .  (20),  PH(x)  can  then 
s  p 

be  approximated  by  rQ(0)  =  tt/4,  yielding 

P 

<Ia>  “  .C  J0(^  y)  e'y  y*y>  2  »  2i 

y  0 

/  2  2  5/3V1 

where  z^  =  ^0.5  k  y  .  Taking  the  upper  limit  as  infinity, 

the  integral  can  then  be  evaluated^)  to  yield 


2  2 
-0/20: 
a 


_  2  2 
8z  9 

a 


D2  e  _p2 / 4qz3 

16qz3 


where  p  *  0z  and 


6  -  (2qz)l/a  »  1.1  C  r1/6z1/a 


n  o 


(27) 


is  the  angle  corresponding  to  the  standard  deviation  of  the  gaussian 
beam  pattern.  Using  Eq.  (27)  and  the  inequality  Eq.  (24)  yields 


e  »  rr 

a  kl 


Hence  the  beam  spreading  in  this  region  is  dominated  by  many  scatter¬ 
ings  from  the  smallest  eddies,  and  for  z  ■  100  km,  6  as  0.2  mrad  for 

a 

typical  turbulence  conditions. 

It  should  be  noted  that  conservation  of  energy  requires 


2tt  J  <I(p,z)>  pdp  -  tt  I 


aperture 


(28) 


which  is  seen  to  be  satisfied  by  Eq.  (26). 

CASE  B:  l  «  p  «  L 
_ ooo 

Substituting  the  spectrum  of  Eq.  (13)  into  Eq.  (9),  expanding  the 
integrand  for  Iq  «  p  «  Lq,  and  integrating  yields 

Mpw(p,z)  -  exp  •[-  1.46  k2  C2  z  p5/3  [1  -  0.80(p/LQ)1 /3  +  0(p/Lo)2]} 

'  °  (29) 

The  second  term  in  the  exponent  proves  to  be  an  important  correc- 

(3) 

tion  to  the  wave  structure  function  given  by  Tatarski,  whose  analy- 

-11/3 

sis  is  based  on  a  spectrum  which  continues  to  increase  as  K  with 
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decreasing  K,  even  for  K  £  .  A  more  complete  discussion  of  this 

point  is  given  in  Ref.  7. 

The  modification  of  Eq.  (29)  for  spherical  waves  can  be  written 
Mg(p,z)  =  exp  {-  0.55  k2  C2  z  p5/3  [1  -  0.71(p/LQ)l/3  ]]■  (30) 

For  Eq.  (30)  to  describe  correctly  the  mutual  coherence  between 
different  elements  of  the  aperture  for  all  p's  of  interest,  it  is  nec¬ 
essary  that  Ms(^0>z)  «  1  and  Ms(Lo>z)  «  1,  which  is  essentially  the 
condition 


z  «  z  «  z.  (31) 

c  i 

Defining  Cl  ■  (0.55  k2  C2  z)3^5  D,  it  follows  that  for  D  Z  «.  , 

n  o 

0  »  1 .  Then  in  Eq.  (20),  r^(x)  «  tt/4,  and  substituting  u  “  fix  re¬ 
sults  in 


<Ib> 


2B2 

Q2 


i: 


S(z/zc) 


3/5 


u  J 


.(n1  u) cxp  " 


5/3 


[‘ 


°'67(-f)1/5“1/3J} du 


(32) 


where,  within  the  range  indicated  by  Eq .  (31),  the  integral  will  be 
insensitive  to  the  precise  value  of  the  upper  limit.  For  D  ~  re¬ 
quiring  M  (D,z)  «  1  leads  to  the  condition 
8 


Z  » 


_ ZJL _ 

1  -  0.71(D/t  )l/a 
o 


(33) 


where  we  have  defined  z. 


(0.4k2  C2  D5/3)_1, 
n 


Substituting  for  0,  Q,  and  a,  one  obtains 


<v 


0.18  D2 

ka/s  ci3/s  zis/6 
n 


(34) 


where  we  have  defined 


1/6 


F(f)  ■  J u  {•  “6/3[l  -  ul/3]} du 


!  u  Jo(£  u)  exp  \ 
b 

[•  “3/3[ 

>  -  °-67(f)  “1/3] 

j-  du 

F| 

(jT 
(z  . 

) 

and 


0.  -  0.69  kl/B  Ce/s  z3/5 
d  n 


(35) 


The  function  F  is  a  slowly  varying  correction  to  the  on-axis 

intensity,  while  the  function  G  determines  the  angular  distribution 

(which  is  "almost"  gaussian) .  F(z/zc)  is  plotted  in  Fig.  8a,  while 

in  Fig.  8b  we  plot  the  result  of  a  mimetic? 1  inversion  of 

G(z/z  ,  9/9,)  =  1/2,  the  latter  giving  the  dependence  of  the  half- 
c  b 

intensity  angle  on  range.  It  follows  from  Fig.  8  that  for 

15  ~  z/z^  ~  104,  one  can  use  the  formula  9i/a  °  0.7  k1^6  C^E  z3^5 

within  an  error  of  w  ±25  percent. 

The  inequality  (31),  together  with  the  condition  z  »  Zp,  implies 
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» 


_1_ 

kD 


J. 

kt 


777  ,  7~Z~  «  0. 


The  scattering  in  the  region  where  Eq.  (34)  applies  is  then  domi¬ 
nated  by  the  presence  of  the  largest  inhomogeneities  and  is  valid  only 
when  6^  is  much  greater  than  the  vacuum  diffraction  angle. 

The  results  of  this  section  are  summarized  in  the  table  on  the 
following  page,  where  the  conditions  required  for  the  use  of  the  ap¬ 
proximate  formulas  are  indicated. 


SUMMARY  OF  MTFs  AND  BEAM  SPREADS 
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V.  GENERAL  DISCUSSION  AND  CONCLUSIONS 


From  Eq.  (18)  one  can  calculate  the  mean  intensity  distribution 

for  an  arbitrary  disturbance  in  an  aperture  if  the  MTF  is  known.  For 

example,  for  a  focused  beam  without  amplitude  perturbations  (i.e., 

2 

-ikr  /2f 

U(r)  =  e  ),  the  mean  intensity  is  readily  shown  to  be  given  by 

kd2 

Eq .  (20),  where  now  in  the  integrand  (only)  0  =  \—  -  j|  and  f  is 

the  focal  length. 

The  MTF  has  previously  been  shown  to  determine  the  limiting  resolu 

tion  obtainable  looking  along  an  atmospheric  path,  as  well  as  the  slg- 

(8) 

nal-to-noise  ratio  of  a  heterodyne  receiver.  '  We  would  like  to  empha 
size  here  that  in  order  to  predict  the  beam  pattern  from  an  arbitrary 
wavefront  in  an  aperture  for  (say)  design  considerations,  it  is  again 
the  MTF  which  should  be  measured  rather  than  specific  beam  patterns. 

To  clarify  this  point,  we  first  observe  that  Eq.  (20)  is  merely 
the  Fourier-Bessel  transform  of  the  quantity  (8/tt)  02Ms(Dx,z)  Tp(x), 
which  can  be  inverted  to  yield 


D2  Tp(p/D)  Mg(p,z)  =  2tt  J  Jq(^-  p)  <I(p,z)>  pdp  (36) 


where  (l(p,z))  is  the  mean  intensity  at  a  distance  p  from  the  axis  at 

range  z.  For  an  arbitrary  disturbance  in  the  aperture,  the  function 

rQ(x)  given  by  Eq.  (19)  would  be  replaced  by  the  appropriate  overlap 
P 

integral  of  Eq.  (18).  Equation  (36)  thus  provides  a  possible  method 

for  determining  "part  of"  the  MTF  from  measurements  of  the  beam  pattern 

However,  because  Tg(p/D)  =  0  for  p  2  D,  inverting  the  intensity 

distribution  can  give  no  information  regarding  M  (p,z)  for  spatial 

s 
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separations  larger  than  the  size  of  the  aperture.  In  particular,  in 
order  to  determine  the  distance  z^  from  a  beam  pattern  measurement,  it 
would  be  necessary  to  have  an  aperture  diameter  greater  than  the  largest 
scale  of  turbtlence,  Lq,  which  might  be  of  the  order  of  meters.  Hence, 
if  one  can  determine  the  spherical  wave  MTF  for  all  spatial  separations 
at  a  given  range,  one  can  infer  the  intensity  distribution  from  an  ar¬ 
bitrary  aperture  distribution  at  that  range,  while  the  reverse  is  not 
true  unless  apertures  greater  than  the  coherence  length  at  that  range 
can  be  constructed.  Even  if  the  beam  pattern  were  measured,  one  would 

first  have  to  construct  M  (p,z)  from  Eq.  (36)  from  the  given  measure- 

s 

ment  to  determine  the  general  response. 

Finally,  although  there  are  many  conditions  under  which  the  spec¬ 
trum  of  Eq.  (13)  is  not  expected  to  apply,  the  intensity  pattern  can 
still  be  determined  from  Eq.  (15)  if  the  MTF  is  known.  For  field  ap¬ 
plications  employing  coherent  optical  devices,  measurements  of  the 
turbulence  spectrum  are  useful  only  to  the  extent  to  which  they  yield 
knowledge  of  the  MTF.  It  is  suggested  that  MTF  measurements  be  made 
and  correlated,  not  with  detailed  spectrum  measurements,  but  with  the 
gross  meteorological  and  topographical  measurements  of  quantities  which 
might  be  observed  to  determine  the  degradation  of  the  higher  spatial 
separations  of  the  MTF,  as  perhaps  wind  speed,  temperature,  and 
properties  of  the  terrain. 
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Appendix 

PROOF  OF  RECIPROCITY 

Let  [U^(r_,  r^) ,  U2(£,  r,2)]  be  field  at  point  r.  due  to  a 
spherical  wave  source  at  (r^,  r2) »  respectively.  Then 

(V2  +  k2n2)  U^r,  r^)  =  -  4Tt6(r  -  r^)  (A-l) 

(V2  +  k2n2)  U2(r,  r.2>  =  -  4rr5(r  -  r2)  (A-2) 

Multiplying  Eq.  (A-l)  by  U^r.,  r^) ,  Eq.  (A-2)  by  U^r,  Tj),  subtracting, 

integrating,  and  using  Green's  theorem,  we  have 


-  -  4tt£u2 (rx ,  fy)  -  U1(r2,  (A-3) 

If  the  surrounding  surface  S  is  taken  as  a  sphere  of  radius  R,  then 
the  surface  integral  can  be  written 

I  Wsr  •  iknU2)  -  “2(sr  •  ik"“i)]  *2dn  <A-4> 

If  we  take  the  spherical  surface  outside  of  the  region  where  n  4  1,  then 

(9) 

the  Sommerf ield  radiation  condition 

liV-  *(Hf  -  iku)  "  0 

which  ensures  that  at  great  distances  from  the  source  the  field  repre¬ 
sents  an  outgoing  wave,  results  in  the  vanishing  of  the  surface  term, 
and  establishes  the  reciprocity  theorem. 
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